Phase behaviour of binary mixtures of diamagnetic colloidal platelets in an external 
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Using fundamental measure density functional theory we investigate paranematic-nematic and 
nematic-nematic phase coexistence in binary mixtures of circular platelets with vanishing thick- 
nesses. An external magnetic field induces uniaxial alignment and acts on the platelets with a 
strength that is taken to scale with the platelet area. At particle diameter ratio A = 1.5 the system 
displays paranematic-nematic coexistence. For A = 2, demixing into two nematic states with differ- 
ent compositions also occurs, between an upper critical point and a paranematic-nematic-nematic 
triple point. Increasing the field strength leads to shrinking of the coexistence regions. At high 
enough field strength a closed loop of immiscibility is induced and phase coexistence vanishes at 
a double critical point above which the system is homogeneously nematic. For A = 2.5, besides 
paranematic-nematic coexistence, there is nematic-nematic coexistence which persists and hence 
does not end in a critical point. The partial orientational order parameters along the binodals vary 
strongly with composition and connect smoothly for each species when closed loops of immiscibility 
are present in the corresponding phase diagram. 



I. INTRODUCTION 

Dispersions of colloidal plateletlike particles, such as 
gibbsiteQ], , montmorillonite [H H| or iron- rich beidel- 
lite d, are susceptible to the influence of magnetic 
fields, since the particles possess nonvanishing diamag- 
netic anisotropy. When a magnetic field is applied to an 
initially isotropic (I) platelet dispersion, the field induces 
orientational order in the system, thus breaking the rota- 
tional symmetry; an orientationally ordered paranematic 
(P) phase results. The paranematic phase has interest- 
ing optical properties, similar to those of the nematic 
(N) phase. When observed through crossed polarisers, 
samples of gibbsite suspensions have been shown to ex- 
hibit field- induced birefringence 7]. Birefringence gradi- 
ents have also been theoretically modelled for a simple 
model system (3 . The effects of a magnetic field on 
montmorillonite platelets were studied in Ref. Q and 
on hematite platelets in Ref. 0. Unlike the gibbsite 
platelets, hematite platelets are ferromagnetic and an I-N 
transition was not observed; rather the authors found a 
clustering effect whereby chains of particles were formed 
due to the platelet-platelet interactions. Such clustering 
has also been observed in simulations 10]. Experimental 
investigations of gibbsite platelets, whereby the suspen- 
sions were exposed to magnetic fields 0, [H| > showed that 
a paranematic phase occurs in these systems. 

The phase behaviour of rods in external aligning 
fields is well-studied, see e.g. KL2lll7| for studies of one- 
component systems. In Ref. |12| the effect of external 
fields on the phase behaviour of rigid rods, freely jointed 



rods and semiflexible rods was analysed with Onsager 
theory [jq : the P-N transition was found to terminate 
at a critical point in all three cases at a high enough 
field strength. In Ref. [l3[ the theories of Landau and 
de Gennes 19] and of Maier and Saupe 20] were used 
to analyse the effects of an applied field in the nematic 
phase. The magnetic field-induced birefringence in solu- 
tions of tobacco mosaic virus (TMV) particles was stud- 
ied in Ref. [3 both experimentally and theoretically (us- 
ing extensions of Onsager theory). In Ref. [l5| the phase 
behaviour of monodispese rods with varying aspect ra- 
tio was studied using the Parsons-Lee scaling [2l|, HH of 
the Onsager functional. It was found that the bifurca- 
tion density decreases with increasing field strength. The 
nematic order of model goethite nanorods in a magnetic 
field was investigated in Ref. 23[ , also using Parsons-Lee 



theory. The goethite rods were modelled as charged sphc- 
rocylinders with a permanent magnetic moment along 
the long axis of the rods. This encourages the rods to 
align parallel to the field at low field strengths. However, 
goethite rods possess a negative diamagnetic susceptibil- 
ity which leads to alignment perpendicular to the field 
at higher field strengths. These competing effects were 
found to yield rich phase diagrams including biaxial ar- 
rangements of the particles. The phase separation in sus- 
pensions of semiflexible /d-virus particles was studied in 
Ref. (l6j, where a P-N phase transition was found. Ef- 



fects of an external field on the isotropic, nematic and 
smectic-A phases of spherocylinders were compared with 
simulations and theory in Ref. with results from 

both approaches being in good agreement. 
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Even when neglecting positionally ordered phases 
(such as columnar and crystal phases), the bulk phase be- 
haviour of binary mixtures of non-spherical colloidal par- 
ticles can be very rich, often including isotropic-isotropic 
(I-I), isotropic-nematic (I-N) and nematic-nematic (N- 
N) phase coexistence, depending on the value of the size 
asymmetry parameter of the two species. The asymme- 
try parameter may quantify the difference in thickness or 
lengths of the two species. An example are binary mix- 
tures of thick and thin hard rods in an external field [24]- 
I2H ] . A general feature of the phase behaviour of binary 
mixtures is a widening of the biphasic region on increas- 
ing the asymmetry parameter. In a certain range of the 
size asymmetry there is typically an I-N-N and/or an I-I- 
N triple point. Coexistence between two nematic states 
may or may not end in a critical point depending on the 
system under study and the value of the asymmetry pa- 
rameter. A well-studied system is the Zwanzig model for 
binary hard platelets, where the particles are restricted 
to occupy three mutually perpendicular directions. This 
was shown to exhibit rich bulk phase diagrams [27l - l29| . 
We recently explored the phase behaviour of binary mix- 
tures of hard platelets with zero thickness and continu- 
ous orientations [301 ] using fundamental measure theory 
(FMT). 

Platelets can be characterized by a diamagnetic sus- 
ceptibility tensor that is diagonal in the platelet frame 
of reference, with components x\\ m the platelet plane 
and X-L normal to it. The diamagnetic anisotropy A\ = 
X\\ — X± 7^ 0) m general, and it may be positive or nega- 
tive depending on the properties of the platelet material. 
For Gibbsite platelets Ax < 0, therefore the platelets 
tend to align with their normals perpendicular to the di- 
rection of the applied field. In order for the platelets to 
align uniaxially in the presence of the field, the samples 
were placed on a central stage and rotated in a horizon- 
tally applied field. In Ref. [8j, FMT was used to study 
the effects of an external field on the phase behaviour 
of monodisperse platelets. It was found that above a 
critical field strength the P-N coexistence ceases to exist 
and the system is homogeneously nematic. In Ref. pjlj . 
van den Pol et. al. have experimentally investigated the 
general phase behaviour of the boardlike goethite col- 
loidal particles (a-FeOOH) in the presence of an external 
magnetic field. The particles were found to align parallel 
to a small magnetic field and perpendicular to a large 
magnetic field; this had already been known since the 
observations of Lemaire et. al. 32] . This effect is due 
to the particles having a permanent magnetic moment 
along their long axis but the magnetic easy axis being 
the short axis. An exciting prospect is that suspensions 
of beidellite platelets, which have a disk-like morphology, 
have recently been shown to undergo an I-N transition 
Q and the nematic phase aligns strongly in the pres- 
ence of an externally applied magnetic or electric field 
0. These platelets possess a positive diamagnetic sus- 
ceptibility and, as such, a simpler experimental setup 
would be required to investigate the P-N transition; the 



platelets are expected to align with their normals parallel 
to the magnetic field. 

Since Rosenfeld's pioneering work [33T - I35I ] there has 
been much interest in the development of FMT for non- 
spherical particles, see e.g. (36[. In the current investi- 
gation we use the FMT of Ref. [30| . which is the mix- 
tures generalization of the theory proposed in Ref. [37]], 
to study binary mixtures of diamagnetic platelets in a 
magnetic field. We consider three different size ratios rep- 
resentative of the different topologies of the bulk phase 
diagram and the full range of external field strengths. 
We investigate how the phase behaviour for each of these 
three size ratios changes on increasing the external field 
strength, which we take to scale with the platelet area 
and to induce uniaxial alignment. 

This paper is organised as follows. In Sec. HI] we out- 
line the density functional theory for the model system. 
The phase diagrams and results for order parameters are 
presented in Sec. [TTT] and we conclude in Sec. HV] 



II. THEORY 

A. Pair Interactions, Model Parameters and 
External Orienting Field 

We consider a binary mixture of hard circular platelets 
with vanishing thickness and continuous positional and 
orientational degrees of freedom. Particles of species 1 
and 2 possess radii R\ and R2, respectively, and we take 
i?2 > Ri- The pair potential uy between two particles 
of species i and j, where i,j = 1,2, models hard core 
exclusion and is hence given by 



Uij(r - y',u),u)') 



00 if particles overlap 
otherwise, 



(1) 



where r and r' are the positions of the particle centres 
and u) and u>' are unit vectors indicating the particle 
orientations (normal to the particle surface). As a control 
parameter that characterises the radial bidispersity we 
use the size ratio 



(2) 



The effect of a magnetic field on the diamagnetic platelets 
is described by an external potential for each species, 



1,2, 



(3) 



where j3 = l/(fcsT), with Ub being the Boltzmann con- 
stant and T absolute temperature; is the angle between 
the platelet orientation u) and the direction of the exter- 
nal field. The strength of the external potential of species 
i is related to the material and field properties via 



Wi 



(4) 
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where B is the magnetic flux density (measured in T) gas of uniaxial rotators is given by 
and Axi 



X 



(0 



is the diamagnetic susceptibility 

anisotropy (with units of JT~ 2 ) of species i, with x^ an d 

x\_ being the susceptibilities perpendicular and parallel 
to the field, respectively 45]. In general both W\ and Wi 
constitute further control parameters. We restrict our- 
selves in the following to special cases and assume that 
Wi scales with the platelet area, i.e. Wi ~ R? . This im- 
plies the relationship Wi = X 2 Wi, and we hence take W\ 
to be our second control parameter, besides the size ratio 
A itself. Scaling with the platelet area is motivated by the 
assumption that the platelets interact with an external 
field in a manner proportional to their mass (neglecting 
any effects of thickness). We could well envisage that 
scaling the strength of the potential e.g. with the radius 
would be another, different yet sensible, choice. We ne- 
glect platelet-platelet interactions due to induced dipoles 
because of their small magnitude, see e.g. the discussion 
in Ref. g|. 

The thermodynamic state is characterised by two di- 
mensionless densities c\ = piR\ and ci = piR\, where 
pi and pi are the number densities of the two species, 
Pi = Ni/V, where iV, is the number of particles of 
species i = 1,2 and V is the system volume. The 
composition (mole fraction) of the (larger) species 2 is 
x = pi/ {pi + pi) and the total dimensionless concentra- 
tion is c = {pi + pi)R\ = c\ + c 2 . 



B. Density Functional Theory 

Density functional theory (DFT) is formulated on the 
level of the one-body density distributions P i(r, u>) of each 
species i. The variational principle [38j j asserts that min- 
imising the grand potential functional f2 yields the true 
equilibrium density profile, 



6Sl([p 1 ,p 2 ],iJ.i,iJ.2,V,T) 



0, i = l,2 



(5) 



where pi is the chemical potential of species i. The grand 
potential functional is given by 

n([pi,pi], mi, H2, V, T) = F id ([p uP2 \,V, T) 

2 

+ F exc ([p 1 ,p 2 ],V,T) + Y^ f dv f dupiiyV^u) 
«=1 J •* 

(6) 

where the spatial integral (over r) is over the system vol- 
ume V and the angular integral (over u>) is over the unit 
sphere. The inter-particle interactions are described by 
the excess (over ideal gas) contribution to the Hclmholtz 
free energy, F cxc ([pi, p 2 ], V, T). We skip the explicit defi- 
nition of the FMT approximation here; this can be found 
in Ref. 30] . The free energy functional for a binary ideal 
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PF id ([p 1 ,p 2 },V,T) = J2 f ^ f du> Pl (r,u>) 
»=i •* •* 

xpnO^r^A?)-!], (7) 

where Aj is the (irrelevant) thermal wavelength of species 
i. 

For bulk fluid states (i.e. with the density distribution 
not depending on r) the orientational distribution func- 
tions (ODFs), ^i(0), i — 1,2, are related to the one-body 
density distributions by P i(r,Lj) = pi^i(9). There is no 
dependence of the ODF on the azimuthal angle 4> since 
the platelets are uniaxial rotators, and we assume that 
only uniaxial states are formed. A powerful feature of 
DFT is that V^ t (r,uj) © appears explicitly in the grand 
potential and therefore enters straightforwardly into the 
minimisation procedure ([5]); see the appendix for the ex- 
plicit form of the corresponding Euler-Lagrange equa- 
tions that we solve numerically. 

The requirements for phase coexistence between two 
phases A and B are the mechanical and chemical equi- 
libria between the two phases and the equality of tem- 
perature in the two coexisting phases (which is trivial 
in hard-body systems). Hence we have the non-trivial 
conditions: the equality of pressure p A — p B and the 
equality of chemical potentials pf = pf , where 2 = 1,2 
again labels the species, and A, B labels the phase. We 
calculate the total Helmholtz free energy F = F^ + F cxc 
numerically by inserting ^i(9) into the free energy func- 
tional. Likewise, the pressure can be obtained numeri- 
cally as p — —F/V + 2j=i Pi9(F/V)/dpi and the chem- 
ical potentials as pt = d(F/V)/dpi. We define a reduced 
pressure p* = fipR\ and reduced chemical potentials 
p* = (3 pi. The equations for phase coexistence are three 
equations for four unknowns (two statepoints each char- 
acterised by two densities) hence regions of two-phase 
coexistence depend parametrically on one free parame- 
ter (which can be chosen arbitrarily, e.g. as the value of 
composition x in one of the phases) and are solved nu- 
merically with a Newton- Raphson procedure 39]. The 
resulting set of solutions yields the binodal. P-N-N triple 
points are located where the P-N and N-N coexistence 
curves cross. In the fieldless case there is, of course, not 

a paranematic phase, but an isotropic phase. 

Pi), 

We characterize orientationally ordered phases (P and 
N) of the binary mixture by two partial order parameters, 
S\ and 52, defined by 



,7T/2 

S l = 4ir d0sin(0)tfi(0)P 2 (cos0), (8) 
Jo 

where P 2 (cos0) = (3 cos 2 6 — 1)/2 is the second Legendre 
polynomial in cos 9. 
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III. RESULTS 

We first review the behaviour of the pure system un- 
der the influence of an aligning field [8]. The inset of 
Fig. HJa) shows the phase diagram for a system com- 
posed of particles of species 1 only. Upon increasing the 
field strength Wx, the coexisting concentrations ci ini- 
tially shift to lower values. The biphasic density gap de- 
creases slightly as the strength of the external potential 
is increased. At approximately Wx = 0.02, the parane- 
matic coexistence concentration starts to increase, while 
the nematic coexistence concentration continues to de- 
crease. The two branches of the binodal meet at a critical 
point at cf tt = 0.42. For Wx > Wf lt = 0.045, there is 
no longer a phase transition and complete destabilisation 
of the P-N transition results in a homogeneous nematic 
phase. For the case of equal sizes of the two components, 
the pure system of species 2 possesses the same phase 
diagram, see the binodal for the case A = 1 in main plot 
of Fig. QJa). However, due to the definition of C2 (recall 
that C2 = P2R1, using the radius of species 1 in order to 
obtain a dimensionless quantity) and the scaling of W% 
with the square of the size ratio (W2 = \ 2 Wi), the phase 
diagram of the pure system of species 2 displays strong 
variation with size ratio A, as shown in Fig.QJa). A shift 
to both smaller values of C2 and of W% occurs upon in- 
creasing the value of A. However, this effect is entirely 
due to the choice of coordinates, which for both pure 
systems are related via C2 = C1/A 3 and W2 = W\/\ 2 . 
Numerical values for the location of the critical point are 
summarised in Tab. HI 



A 


A 2 


A 3 


C C 2 rit 


WT rit 


1 


1 


1 


0.42 


0.045 


1.5 


2.25 


3.375 


0.124 


0.019 


2 


4 


8 


0.053 


0.011 


2.5 


6.25 


15.625 


0.027 


0.0072 



TABLE I: Scaling of the location of the critical point with 
size ratio: Critical concentration c§ nt = A _3 c crlt and critical 
field strength Wi crit = A" 2 W crit for a range of size ratios A, 
where c crlt and W crlt are the critical concentration and field 
strength in the pure system (without species index). 

The variation of the order parameter S for monodis- 
pcrse platelets with field strength is displayed in 
Fig. [1(b). I n the field- free case, W\ = 0, the nematic 
phase at coexistence possesses an unusually small order 
parameter, see e.g. the discussion in Ref. [4(|. For all 
size ratios, as Wx is increased, the coexistence value of 
S in the paranematic phase increases monotonically, and 
the value of S along the nematic branch of the binodal 
decreases with increasing field strength. This is consis- 
tent with the fact that the coexistence density decreases 
as the field strength increases, overcompensating for the 
ordering effect caused by the applied field. At the criti- 
cal point the nematic order parameter takes on the value 



52 = 0.27. For increasing values of A = 1, 1.5, 2 and 2.5, 
the critical point shifts to smaller values of field strengths. 

We next consider the binary mixture in the external 
field and hence explore the full range of compositions, 
< x < 1. In Fig. [2] the phase diagram for A = 1.5 is 
shown. We consider a range of external field strengths 
up to Wx =0.15 (which corresponds to W2 — X 2 Wx = 
0.3375). For the fieldless case, Wx = 0, there is FN phase 
coexistence over the entire range of compositions x. We 
display this phase diagram (and subsequent ones) both 
in the (cx,C2) representation [Fig. 12(a)] as well as in the 
(x,p*) representation [Fig. 11(b)] . Tie-lines are omitted 
for clarity; in the (01,02) representation these connect 
the lower branch of the binodal to the upper branch in 
such a way that the isotropic (or paranematic) phase is 
rich in (the smaller) species 1 and the nematic phase is 
rich in (the larger) species 2. In the (x,p*) representation 
[Fig. E(b)] the tie lines are (trivially) horizontal due to 
the condition of equal pressures in the coexisting phases. 
For Wx — 0.01, the binodal still connects to the axes 
(which correspond to the pure systems). Recall that the 
P-N transition still occurs in the pure systems at this field 
strength, cf. Fig.QJa). However, the isotropic phase has 
now become a weakly-ordered paranematic phase. Hence 
there is P-N phase coexistence over the entire range of 
compositions. The isotropic phase has been replaced by 
a paranematic phase, because the order parameter along 
the lower branch of the binodal is non-zero, see Fig. [2(c) 
and (d), where the partial nematic order parameters are 
shown for species 1 and 2, respectively. On increasing 
Wx to 0.03, the upper and lower branches of the binodal 
still persist to the pure system of smaller platelets, con- 
sistent with the findings of Ref.Q. However, the binodal 
does not touch the C2-axis, indicating that there is no 
longer a phase transition in the pure system of species 
2 (we found the critical field strength for the monodis- 
perse system at A = 1.5 to be 0.019, which is less than 
0.03, Fig. {I]})). Hence the two branches of the binodal 
connect at a (lower, in pressure) critical point. There- 
fore the state of the system changes continuously from 
paranematic to nematic for compositions greater than 
about 0.7 by increasing the pressure. For compositions 
less than this value, increasing the pressure from below 
the lower branch of the binodal to the upper branch of the 
binodal, the system, as before, passes through a bipha- 
sic region. For Wx = 0.05 the departure of the binodal 
from ci = (and x = 0) occurs as is consistent with 
the critical field strength being W^ rlt = 0.045 in the pure 
system. The result is a phase diagram in which the two 
branches of the binodal have joined to form a closed loop 
of immiscibility. There is a larger range of compositions 
towards the x = 1 side of the phase diagram (approx- 
imately x > 0.45) than towards the x — side of the 
phase diagram (approximately x < 0.02), where an in- 
crease in pressure leads to a continuous change in state 
from a paranematic state to a nematic state. The order 
parameter of species 1, measured along the binodal varies 
with composition [Fig. (J3b)] such that for Wx = 0, 0.01 
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FIG. 1: Paranematic-nematic phase diagram of the one- 
component system(s). (a) Behaviour of the pure system of 
species 2 (i.e. x = 1). The variation of the paranematic (P) 
and nematic (N) coexistence concentrations, C2 (horizontal 
axis) with the strength of the aligning field Wi (vertical axis) 
is shown for A = 1, 1.5,2 and 2.5 (from right to left, as in- 
dicated). The inset shows the phase diagram for the pure 
system of species 1 (i.e. x = 0). This is equivalent to the case 
A = 1 in the main plot when identifying the horizontal axes. 
The critical point is depicted as a filled circle, (b) Variation 
of the orientational order parameter S — S2 along the parane- 
matic and nematic branches of the binodal (horizontal axis) 
of the pure system of species 2, with increasing field strength 
W\ (vertical axis) for the same size ratios as in (a). Critical 
points are depicted as filled circles. 



FIG. 2: Phase behaviour of binary platelet mixtures with 
size ratio A = 1.5 and external potential strength Wi = 
0,0.01,0.03,0.05,0.07,0.1,0.12 and 0.15 (with W 2 = \ 2 Wi) 
from outside to inside. Shown are phase diagrams in (a) the 
(ci,C2) representation and (b) in the (x,p*) representation. 
The partial order parameters Si along the binodal are shown 
in (c) for species 1 and in (d) for species 2. 



and 0.03 the paranematic and nematic branches of the 
binodal do not connect on the low-composition side of 
the order parameter graph since these values of W\ are 
less than Wf x% = 0.045. For 0.05 < Wi < 0.15 the two 
branches of the binodal connect on the low-composition 
side of the graph. For W\ = and 0.01 the two branches 
of the binodal do not connect on the high-composition 
side of the graph since these values of W\ are smaller 
than VF 2 crit = X 2 Wf lX = 0.019. For W x > 0.05 the par- 
tial order parameters measured along the binodal form 
closed loops. These islands become smaller with increas- 
ing field strength and eventually coalesce to a point when 
the double critical point is reached. The partial order 
parameters of species 2 [Fig. ^\)] follow a similar pat- 
tern except that the order at a given statepoint is higher 
than that for species 1, as one could expect, given that 
species 2 is of the larger size. 

In Fig. [3] we show results for A = 2. Increasing the 
size ratio to this value leads to an increase of the size of 
the I-N biphasic region [30(. The fieldless case possesses 
a reentrant phenomenon whereby the system undergoes 
the following change of state when increasing the pres- 
sure at fixed mole fraction at around x — 0.7 starting 
in the isotropic region: I— >I-N— >N^I-N— >N. In addition, 
there is N-N coexistence between a nematic phase rich 
in species 1 (Ni) and a nematic phase rich in species 2 
(N2) ending in an upper critical point and an I- N-N triple 
point. Applying a small field strength of W\ — 0.02, the 



FIG. 3: Same as Fig. [2] except for size ratio A = 2 and for 
external field strengths Wi = 0, 0.02, 0.6, 1 and 1.2. 



binodal no longer reaches the pure system of species 2, 
as Wi = 0.02 > Wf 1% = 0.011, which is the critical field 
strength for the monodisperse system at A = 2. An effect 
of this is that the reentrant part of the phase diagram 
alters: the range of compositions for which the system 
undergoes P— >P-N— >N— kP-N— >N at just over x = 0.5 is 
much smaller than in the fieldless case. The binodal ends 
in a tail-like feature at just under x — 0.6. Aside from 
these differences, the rest of the phase boundaries fol- 
low closely those of the fieldless case, though remaining 
slightly inside those of the latter. There is N-N coexis- 
tence ending in an upper critical point. The triple point 
is retained as a P-N-N line in the (x,p*) representation 
and a triangle in the (01,02) representation although we 
do not show these features in the plots for clarity. Upon 
increasing the field, triple phase coexistence vanishes, i.e. 
the triple point collapses onto two-phase coexistence. We 
have not calculated the precise value of the external field 
where this happens. We expect this value to be different 
from the values where the binodal detaches from either of 
the density axes (i.e. differ from the critical field strengths 
in the pure systems). Applying a field strength W\ = 0.6, 
which is much greater than the critical field strength for 
the monodisperse case, W" lt = 0.045, leads to a closed 
loop of immiscibility. The nematic phase rich in small 
platelets (N2) and the paranematic phase have merged. 

In Fig. 2] we present results for A = 2.5. The field- 
less case displays TN coexistence, an I-N-N triple point 
and coexistence between two nematic states, which does 
not end in a critical point. Applying just a small field 
W\ = 0.01 has a considerable effect on the phase be- 
haviour: the pure system of species 2 looses the P-N 
transition as W\ = 0.01 > Wf 1 * = 0.0069. However, 



FIG. 4: Same as Fig. [2] except for A = 2.5 and for external 
field strengths Wi = 0,0.01,1,2.5,5 and 7. In (c) and (d) 
we only present partial order parameters for Wi = 0.01, 1, 2.5 
and 5 for clarity. The small kinks at low values of Si in (c) 
and (d) are due to numerical artifacts. 



the transition persists in the pure system of species 1. 
Hence, there is still P-N coexistence and indeed a tail 
between about x = 0.5 and x = 0.6 where reentrant 
behaviour occurs. For W\ — 1, there is no P-N tran- 
sition for either of the pure limits. There is, however, 
a large immiscibitity gap between two distinct nematic 
phases, N\ and N 2 . Increasing the field strength raises 
the phase coexistence to higher pressures and narrows 
the phase coexistence region. However, large steps in 
field strength are required to have a significant effect on 
the system. W\ = 7 (which corresponds to W% = 49) 
is approximately 45 times stronger than the field rcqircd 
to homogenise the system at A = 1.5 and yet even with 
such a high field strength, there persists a wide coexis- 
tence region. The order parameter curves approach very 
high values (close to 1) for large values of W\. This forms 
a limit to the densities at which we may probe at A = 2.5: 
for very high order parameters it is numerically difficult 
to determine the ODFs, even with very fine #-grids. 



IV. CONCLUSIONS AND OUTLOOK 

We have investigated the effects of an external align- 
ing (magnetic) field on the phase behaviour of binary 
mixtures of circular hard core platelets with zero thick- 
ness. Using the FMT of Refs. 0, we have traced 
paranematic and nematic phase boundaries and have ex- 
amined the partial nematic order parameters at coexis- 
tence. Three different representative values for the ra- 
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dial bidispersity (A = 1.5, 2, 2.5) have been studied. The 
topologies for each of these values are different from each 
other in the fieldless case. For the smallest size ratio con- 
sidered, A = 1.5, the fieldless case shows only I-N phase 
coexistence over the entire range of compositions. For 
A = 2, besides FN coexistence there is also N-N coex- 
istence ending in an upper critical point and an FN-N 
triple point occurs. For A = 2.5 the N-N phase coexis- 
tence does no longer end in a critical point (at least up to 
the densities we considered). Applying the external field 
induces paranematic order in the low-density regions of 
the phase diagrams (which are isotropic in the absence of 
a field) . Increasing the field strength leads in all cases to 
a narrowing of the biphasic P-N region. The P-N transi- 
tion further destabilises upon increasing the external field 
strength. The system becomes more strongly ordered, 
such that for A = 1.5 and 2 phase coexistence disappears 
(at a double critical point) and the system is in a single- 
phase nematic state for all statepoints. The field strength 
required to complete the homogenisation increases for in- 
creasing size ratio. In contrast, for A = 2.5, the coexisting 
nematic states become so well-ordered that the system 
does not become homogeneously nematic up to the field 
strengths we have applied. 

Results from computer simulation studies for this (or a 
similar) model mixture are highly desirable, as are exper- 
imental studies. Colloidal platelets are often significantly 
polydisperse in both radius and thickness, see for exam- 
ple Ref. [4l|, so effects due to polydispersity will play a 
role in experimental systems, which are not accounted 



for in the present theory. Recently, the P-N interface in 
suspensions of boardlike goethite particles has been in- 
vestigated experimentally [42| . Anticipating that similar 
studies could be made in systems of colloidal platelets, a 
further exciting avenue would be to investigate the prop- 
erties of the P-N interface using FMT, which has already 
been shown to compare well with experimental and sim- 
ulation results for the I-N interface [43]. It would be in- 
teresting to consider in theoretical work the effects that 
are induced by finite thickness of the platelets. In the 
present study we restricted ourselves to mixtures with 
moderate size asymmetry, as we expect the theory to de- 
scribe these accurately. Investigating highly asymmetric 
mixtures, possibly based on the depletion picture, is an 
interesting issue for future work. 



Appendix A: Self-Consistency Equations for the 
Orientational Distribution Functions 

We give a summary of the equations that are necessary 
to find the ODFs at a given composition, x and conce- 
tration, c. The excess free energy from FMT is the sum 
of the right hand sides of Eqs. (13) and (27) of Ref. [3(j. 
Inserting this, together with the ideal free energy (JTJ) and 
the external potential ([3]) , into the grand potential func- 
tional ([6]), and employing the minimisation principle ([5]) 
leads to two coupled Euler-Fagrange equations for the 
ODFs: 



tfi(0) = ^exp 



,7r/2 1 

8ttc / d9' sin 9'K(9, 9') [(1 -x)Vi {9') + -x(\ 2 + \)^ 2 (9')} 
Jo 2 

tt/2 r ir/2 



- 327TC 2 / d9'sm9' d9" sm9"L(9, 9', 9") 
Jo Jo 

[(1 - a:) 2 * + 2a?(l - x)\ 2 V 1 {8')t> 2 {9") + x 2 X i ^ 2 (9')^ 2 (9")} + W x sin 2 i 



(Al) 



*a(0) = ^exp 

^2 



prr/2 

32ttc 2 / d9 
Jo 



,tt/2 , 

8ttc / d9' sm9'K(9,9')[x\ 3 $ 2 (9') + -(1 - x)(\ 2 +A)#i(0 / ) 
Jo 2 

prr/2 

'sm8' / d6" sin 9"L(9, 9', 9") 
Jo 



[x 2 \ 6 ^ 2 (9')^ 2 (9") + 2a?(l - x)X 4 ^ 1 (9')^ 2 (9") + (1 - a;) 2 A 2 *i(0')#i(0")] + W 2 sin 2 9 

I 



where the constants Z\ and Z 2 are such that the nor- The integral kernel K(9,9') is 

malisation J du:^i(u)) = 1, for i = 1,2. The numerical 27r 27r 

procedure is the same as that described in Ref. [3(| , which K{6, 9') — dip sin 7 = / d<fr-\/l — {u> ■ u) 1 ) 2 

is an extension of the procedure introduced in Ref. [44| . jo Jo 

2 ir 



d4>^/\ - (cos 9 COS 9' + 



(A2) 



cos cos tr + sin v sin fr cos < 



(A3) 
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where <fi is the difference between the azimuthal angles of 
the two platelets and the kernel L(9,6',9") is 



L(9,9',9") = 



2-K 



JO 
2-k n 2n 



\u ■ («' X u") 



cos 9' sin <f>" sin 0") + cos 6>(cos 4> sin 9' sin 4>" sin 9" 
sin 4> sin 9' cos 0" sin 0") | . 



(A4) 



The solutions of Eqs. flS) and ([A2]) . * x (0) and * 2 (0), 
are then inserted into the equations p A = p B and /if — 
/if, i = 1, 2 in order to find the coexisting states. Triple 
points are located where the P-N binodals and the N-N 
binodals intersect. 
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